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Suppression of Tunneling in a Superconducting Persistent-Current Qubit 
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We consider a superconducting persistent-current qubit consisting of a three-junction supercon- 
ducting loop in an applied magnetic field. We show that by choosing the field, Josephson couplings, 
and offset charges suitably, we can perfectly suppress the tunneling between the two oppositely 
directed states of circulating current, leading to a vanishing of the splitting between the two qubit 
states. This suppression arises from interference between tunneling along different paths, and is 
analogous to that predicted previously for magnetic particles with half integer spin. 

PACS numbers: 



Surprising physical effects can be produced by 
quantum-mechanical interference between particles mov- 
ing from one site to another along different paths 0,0. 
Examples include the Aharonov-Bohm effect Q and the 
Aharonov-Casher effect 4, 5]. Many authors have consid- 
ered this Aharonov-Casher effect for Josephson-junction 
arrays and devices [E H, IE El- Starting from this topolog- 
ical effect Loss et al.[l(| predicted suppression of tun- 
neling due to interference of different tunneling paths 
for magnetic particles with half integer spin, and also 
oscillations in tunnel splitting with applied magnetic 
field. These oscillations were confirmed experimentally 
by Wernsdorfer and Sessoli [ill ]. 

Effects related to those analyzed in Ref. [l(| have also 
been studied in systems of Josephson junctions [1, [H, [H|] • 
For example, a three-junction loop has been studied as 
a possible phase qubit [l2l. [l3j. In its original design [T^|. 
the junctions were deliberately made asymmetric to avoid 
interference of different tunneling paths, and to protect 
the qubit from random charge fluctuations. 

In this note, we analyze the same three-junction loop 
but in a different regime, namely, one in which the asym- 
metry is chosen to observe and control the interference of 
tunneling paths. We show that for this chosen asymme- 
try, the interference is easily detected because the tun- 
nel splitting vanishes perfectly for certain special induced 
gate charges. The suppression of the tunneling splitting 
is closely analogous to that discussed in Ref. [l(|. Al- 
though this qubit is potentially subject to random charge 
noise, it may be possible to minimize this noise by ap- 
propriately tuning the gate voltage [l4j. 

The circuit diagram for the three-junction qubit of Ref. 
[13] is shown in Fig. [Q The i th junction (i = 1, 2, 3) 
has capacitance C\ and Josephson coupling energy Eji. 
An external flux $ = f^a, where $o = h/(2e), is ap- 
plied through the three-junction loop, which is assumed 
to have negligible self-inductance. The properties of the 
qubit can be manipulated by controlling <&, and also the 
two external voltages, Va and Vb, which are applied to 
the circuit through two gate capacitors C 9 a and C g s- 
The voltages across these capacitors are V q a = Va — V\ 
and V gB =Vb -V 2 . 



We assume, following Ref. that junctions 1 and 2 
have equal Josephson energies Ej and equal capacitances 
C, while junction 3 has Josephson energy aEj and ca- 
pacitance aC. We also assume that C 9i a = Cg,B = jC. 
Because of the fluxoid quantization, the three gauge- 
invariant phase differences (pi, <f>2, and 03 satisfy the 
constraint <f>\ — 02 + 03 = —2irf. We choose 0i and 

2 to be the independent dynamical variables, and write 

03 = 27r/ + 4>i - 4> 2 - 

With these choices, the circuit Lagrangian C can be 



written as C = T - U - (§a)0 • C g ■ V g = £ Q + C W z, 
where Cq = T — U and C\yz is the remaining term, 
which we call the Wess-Zumino term. The kinetic en- 
ergy term T represents the electrical energy stored in 

all the capacitors of the system; it can be written as 
—T — 

T = \{2^) 2 <t> ■ C ■ <t> The Josephson coupling energy 
U = Ej[2 + a — cos 0i — cos ^2 — acos(27rf+ <j>\ — $2)]- 



■ C g ■ Vg. 



Finally, £ wz = 

In the above expression for C, the quantity C repre- 
sents the capacitance matrix for the junctions, with ma- 
trix elements Cn = C22 — C(l + a + 7), C12 = C21 — 

—aC. Also, </> represents the (1 x 2) column vector with 
elements (j>i, 4>2] V g is the 1x2 column vector with el- 
ements Va and Vb', and the gate capacitance matrix C g 
is a 2 x 2 diagonal matrix with both diagonal elements 
equal to 7C. All the above formalism is identical to that 
in [T2I ]. However, we consider a different set of param- 
eters: a > 1 (rather than a < 1) and / ~ 1/2. This 
difference has a qualitative effect on the qubit behaviour. 

Fig. [2] shows a contour plot of the potential energy 
^(01,02) for ct — 1.3 and / = 1/2, represented in the 
repeated cell scheme. The origin ((f>i = 0, 2 = 0) is 
shown at the center of the plot. The potential energy is 
periodic in both tpi and <f> 2 with period 2-7T. With our 
choice of parameters, this potential energy has two in- 
equivalent states of the same minimum energy, indicated 
by boxes with horizontal and vertical lines within each 
unit cell. If we choose the unit cell to be — n < <j>x < 7r, 
— 7r < <p2 < 7T, then the two inequivalent minima are 
in the upper left and lower right quadrants of the cell. 
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These two minima occur at (0* + 2to7t, — 0* + 2nir) and 
(— 0* + 2m7r, 0* + 2nn) where m and n are positive or 
negative integers and 0* = cos -1 These two states 
are degenerate, but inequivalent. Physically, they corre- 
spond to states with clockwise and counterclockwise loop 
currents. When a > 1, the lowest-barrier tunneling paths 
between these states are different than in the a < 1 case 
considered in Ref. [12j |. 

To see this, suppose we start from the state (0* , —(f)*), 
and suppose that this state represents a clockwise- 
circulating loop current, corresponding to boxes contain- 
ing vertical lines in the lower right-hand corner of the cen- 
tral unit cell. There are three plausible tunneling direc- 
tions to reach a neighboring state with counterclockwise- 
circulating currents, leading to states at (2n — 0*,0*), 
{-<fi*, -27T + 0*), and (-0*, 0*). If a > 1, one can show 
numerically that the potential barrier is smaller for the 
two located at (2ir — 0*,0*) and (— 0*, — (2tt — 0*)) than 
that for tunneling to (— (f)*,(j>*). By contrast, if a < 1, the 
tunneling barrier is smaller for the third path than for the 
other two. Because there are two possible lowest-barrier 
paths when a > 1, there is an interference effect in this 
case which is absent when a < 1. Furthermore, the dif- 
ference in barrier heights (between the two equal-barrier 
tunneling paths and the third, higher-barrier path) in- 
creases with increasing a, provided a > 1. Thus, we can 
easily choose a so that the system tunnels only through 
these barriers. This tunneling corresponds to the paths 
near the heavy line in Fig. [2] We will show that, for 
a > 1, there exist certain values Qa and Qb of the 
stored charge, for which the tunneling along these two 
equal barrier paths exactly cancels out. 

In the absence of tunneling, the system has two degen- 
erate minimum-energy quantum states when / = 1/2, 
one with counterclockwise (0) and the other with clock- 
wise (0) current. In the presence of tunneling, these 
two states are connected by a tunneling matrix element 
w, which breaks the degeneracy. 

The transition amplitude (V) from a state to state 
can be calculated using the imaginary time coherent 
state path integral method. Symbolically, at temperature 
T = 0, V = !i%];tl^ DQe- 1 ^ , where DSl repre- 
sents an integral over all paths in imaginary time starting 
from the clockwise state at (0i(O), 02(0)) = (0*,— 0*) 
at r = and ending at the counterclockwise state at 
(0i(oo),0 2 (oc) = (2tt - 0*,0*) or (-0*, 2?r - 0*) at 
t = oo. S r j >ll j >2 represents the action calculated along each 

of the paths. In turn, S 0102 = J^^^io)^ dr i^o+C wz ), 
where the integral is over imaginary times r (such that 
t = it), and Cq + C wz is the Lagrangian but with each 
time t replaced by it. At T = 0, the integrals start at 
t = and run to r = oo. 

The key point is that, for a > 1, there are two classes 
of paths going from the point to the point in 
phase space. One of these is generally in the "north- 
east" (NE) direction and the other in the "southwest" 
(SW) direction; the paths run in generally opposite di- 



rections in the vicinity of the heavy black line in Fig. 
2. The two end-points of the paths in the NE direction 
are = (0*, -0*), and Q = (2vr - 0*,0*), while for 
those in the SW direction they are = (0*,— 0*) and 

= (—</>*, — (2ir — (j>*)). Let us consider one particular 
path in the NE direction, and denote it by (0i(t), 02 (r)). 
This path runs from (0*, -0*) to (2tt-0*, 0*). Then the 
path — 02 (t), — 0i (t)) also starts from (0*, —0*) but runs 
generally in the SW direction to (-0*, — (2tt— 0*)). Thus, 
for every path in the NE direction, we can define a cor- 
responding path in the SW direction by this procedure. 

We now show that the contributions of these two paths 
to the path integral exactly cancel out for special values 
of Qa and Qb- We first consider the contributions of 
IA and T to the path integral. At any point along a 
NE path, the potential energy U{4>\, 02) > is given above. 
Along any point along the corresponding SW path, the 
corresponding potential energy is given by U{— 02, — 0i), 
which from eq. (3) equals W(0i, 02). Thus, the contribu- 
tion of U to S is exactly the same for corresponding paths 
in the NE and SW directions. Similarly, the contribution 
of T to S is the same for corresponding paths in the NE 
and SW directions (because T is quadratic in the deriva- 
tives 0i and 02 > and because the diagonal elements of C 
are equal). Since S appears in the exponential, the ex- 
ponential exp(— S/H) terms give the same multiplicative 
contribution to V for each of the two paths. 

For Cwz, we have C W z = -i (§£) iC(Va4>i + Vb0 2 ), 
where 0i and 02 are derivatives with respect to r. 
Since Cwz is a total time derivative, its contribution 
to S(f, lt f, 2 depends only on the initial and final values 

1 and / of the phases (<f>i and 02). From eq. (4), 

this contribution is — f. oIt^C^a^i + Vb0 2 )) = 

-i(thc[v A (Mf)-Mi)) + v B (Mf)~Mm- Thus, 

for any path taking the state (0*,— 0*) in the NE di- 
rection to (27r — 0*,0*), Cwz gives a contribution to 
S equal to S W E Z = -i (fa) 7 C [2V A (ir ~ 0*) + 2V S 0*]. 
Similarly, a path taking the state (0*,— 0*) in the SW di- 
rection to (— 0*,— (2tt — 0*)) gives a contribution = 
i(fa) 7 C[2Vk0*+2F B (7r-0*)]. 

We now write §%lCV A = h%± and ^CV B = h%f, 
where Qa and Qb represent the charge stored on the 
gate capacitors, and e represents the electronic charge. 
Then the sum of the contributions of 5m| and S^f^ to V 
is Vwz = "P (e = [9 ' 4(7r ~^* )+Q£! '**] 4- e ~e[ Q ^* +<3B(7r ^*)]), 
where Pq is a constant term which is the same for the two 
paths. It is always possible to define a number n such 
that Qb = ne — Qa- In terms of n, we can rewrite Vwz 
as Vwz = P o exp[iQ A (7r-20*) + m0*] [l + e'™]. 
Hence, Vwz vanishes whenever Qa + Qb = ne, where 
n is an odd integer. Since the contributions of U and 
T to S are the same for corresponding paths in NE and 
SW directions, the total V still vanishes when this con- 
dition is satisfied, even including the contributions of U 
and T to S. Our calculation is analogous to that of Loss 
et a/.[To| for a magnetic tunneling problem. 

Thus, the paths taking the state (0*,— 0*) to (2ir — 
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0*,0*) and to (— 0*,— (2n — 0*)) interfere completely de- 
structively whenever the stored charges on the gate ca- 
pacitors sum to an odd multiple of e. This destructive 
interference is not restricted to straight line paths, be- 
cause, as we have shown, for any general path in the 
NE direction, there exists a path in the SW direction 
which interferes destructively with it. It can be shown 
that similar destructive interference occurs for higher or- 
der paths, such as the next order paths which take the 
state (0*,-0*) to (4tt - 0*, 2tt + <jf) and (-(2tt + 0*), 
— (47r — 4>*)). This implies that the two persistent cur- 
rent (clockwise and counterclockwise) states remain de- 
generate for these values of the stored charges, provided 
that / = 1/2. 

The same cancellation can be demonstrated using the 
tight-binding formulation of Ref. [T^. The classical 
Hamiltonian TL corresponding to the above Lagrangian is 

TL = p T ■ <j) — C, where the canonical momenta pi = i 

= 1, 2. This procedure gives' H = \ (j^F? • CT 1 • 

P + U((pi, 02), where P = p + f^C g • V g , and U is given 
above. 

The energy eigenstates of TL satisfy the time- 
independent Schrodinger equation TLif>(<j)i, 2 ) — 
Ety(0i,02), where p t = -ih^/dfa) (i = 1, 2). The 
boundary conditions, obtained from the require- 
ment that the wave function be single- valued, are 
if>(<j)i + 2to7t, 4>2 + 2mr) — ipi&ii 02), where m and n are 
integers. 

The charge periodicity discussed above is due to the 
relation between the components of P and p. Specif- 
ically, using the fact that the matrix C g is diagonal, 

we write P x = Pl + ^C g , A V g , A = -ih (jjfc + i% ) , 

P2=V2 + ^C g<B V g>B = -ih + ifia ) . on using the 
operator forms of p\ and P2 and using Q A = C gA V g ^A 
and Qb = C g ^B,V g ,B- We now define x(0i> 02) = 
exp(zk • (f>)ijj(<pi, 4>2), where the vector — 0101 + 0202, 
k = ^01 + ^f-02, and 0i and 02 are unit vectors in the 
0i and 02 directions in Fig. 2. Then the boundary con- 
ditions on if> imply that x(0i>02) is a Bloch function, i. 
e., that x(0i + 27rm, 02 + 27rn) = exp(ik- R mn )x(0i, 02), 
where R = 27tto0i + 27rn02 is a lattice vector. Also, in 
terms of x, the two-variable Schrodinger equation takes 
the form (T +U)xu = i?(k)xk- Since Xk is a Bloch func- 
tion and U is periodic, the eigenvalue E(k.) is periodic in 
k and hence, in the charges Qa and Qb. 

This Schrodinger equation can be solved within a tight- 
binding approximation [12| to calculate the tunnel split- 
ting, which we previously calculated using a path inte- 
gral approach. Let us consider two localized "atomic" 
orbitals, w(0i,02) and w(0i,02), which represent the 
ground state wavefunctions in each of the two minima 
of the potential W(0i,02) within the central unit cell. 
Then, at Bloch vector k ,the tight binding wavefunction 
Xk(0i,02) = c k , ll w(0i,02) + Ck, 1 >w(0i,02), where c k ,u and 
Ck >v are defined by [H uu (k) - B(k)]c k ,« + H uv (k)ck <v = 0, 



-fftm(k)c k ,u + [H vv (k) - E(k)]c k . v = 0. In actuality, H uu 
and H vv are independent of k. 

When the applied field is such that / = 1/2, u and v 
are exactly degenerate, with energy which we denote eo- 
In this case, H uu = H vv = €q. To obtain the other two 
elements, we denote by t\ the tunneling matrix element 
between these two minima in the same unit cell, and 
t2 between the state in the "southeast" corner of that 
cell and either of the two adjacent minima lying along 
the heavy line in Fig. 2. Then H uv = H* u = -t x - 
t 2 [e 4k ' Rl + e 4k - R2 ] , where R x = 2tt0i and R 2 = -2tt0 2 
are the Bravais lattice vectors from the central unit cell 
(denoted by a white square in Fig. 2) to the two adjacent 
cells along the heavy line in Fig. 2. 

An estimate of the t^s can be obtained using the WKB 
method, by calculating the action Si between the two 
minima and writing tj ss {hiOi/2iT)e^ s ^ h . Here u>i is 
the attempt frequency for the phase "particle" to escape 
from the potential well. Following the approach of Ref. 
P3, we find that for a = 1.3 and E, ; /E c ~ 100, where 
E c = e 2 /(2C), the ratio t 1 /t 2 ~ 10" 4 . In fact, provided 
Ej/Ec <C 1, we can always choose an a > 1 such that 
the effect of t\ is very small. 

Neglecting t\ and using the above values of k, 
Ri and R2, we obtain H uv = — 2t2 exp[z7r(<5A — 
Q B )/{2e)] cos[tt(Q a + Qs)/(2e)]. The eigenvalues of H 
are then £^ = 60^ |^uw|, with corresponding normalized 
eigenvectors (u ± v)/y/2. The result for the eigenvalues 
shows that, when the offset charges satisfy Qa+Qb = ne, 
with n an odd integer, the levels become degenerate. This 
is exactly the result we found by our path interference 
analysis. Note that the energy splitting depends only on 
the sum Qa + Qb, not on the difference Qa — Qb- 

As / deviates slightly from 1/2, the potential U 
changes such that the two minimum states of the two 
wells in a unit cell become unequal in energy, and the 
barrier heights also change. If we define the zero of en- 
ergy as the average of the two lowest energy states at 
/=l/2, the elements of H become H uu = —H vv = F; 
H U v = H* u = -t. Here F ~ (dH uu /df){f - 1/2) is 
the change in the diagonal matrix element of H with 
a small change in flux. Also, if we write Qa + Qb = 
(2n + l)e + 5Q, and we assume \5Q/e\ <C 1, we find 
i~ -2t 2 exp[in(Q A -Q B )/(2e)](-l) n Tr6Q/{2e). Thus, 
to first order in Sf and SQ, F and t are controlled by 
two different parameters: / — 1/2 and SQ/e. The corre- 
sponding eigenvalues of H are E T = tV f2 + and 
depend on / — 1/2 (through F) and Qa + Qb (through 
t) , but not on Qa — Qb ■ (The eigenvectors do depend on 
Q A — Qb-) By manipulating these two control param- 
eters independently, one could, in principle, adjust the 
splitting of this two-level system. 

Fig. [3] shows a contour plot of the quantity 
2^/F 2 + \t\ 2 , which represents the energy difference be- 
tween the two lowest eigenvalues of H, as a function of 
the quantities / and (Qa + Qb)/ c - ln constructing this 
plot, we assume the following parameters: Ej/Ec — 
80, 7 = 0.02, a — 1.3, and an attempt frequency 
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Huj/(2tt) = 0.193Ej. Except for a, all these quanti- 
ties have the same values as in Ref. [l2j]. There are 
eight contour curves visible, equally spaced between 
and a maximum value of 0.3266-Ej. These are cal- 
culated using the above parameters and the relation 
[dH uu /df] f=1/2 = 4naE j sin(2Trf + 2<j>*). The split- 
ting vanishes when (Qa + Qs)/e is an odd integer and 
/ = l/2. " " 

The parameters F and |f| should be controllable ex- 
perimentally. F can be finely adjusted by changing /, 
the ma gne tic flux through the loop. For the parameters 
of Ref. 12], \t\ should also be controllable. Taking Ej = 
800/ueV, we have Ec — 10/ieV, corresponding to a junc- 
tion capacitance C — 10/F, and hence gate capacitances 
C g = 0.16/F. With this value for C g , (Q A + Qs)/e = 1 
corresponds to Va + Vb — C 9 (Qa + Qb) — lmV, a value 
which should be tunable to a small tolerance. In the dif- 
ferent regime of small junctions (Ej ~ Ec), the periodic- 
ity of e nerg ies with offset charges has been observed, e. g., 
in Ref. 14]. In that work, a computer-controlled method 
was used to accurately compensate for the random offset 
charges. For operation of the present system as a qubit 



in the regime with a > 1, one would need a temper- 
ature T low enough to avoid creation of single-electron 
excitations, or exciting the system above its two lowest 
levels. Using the parameters of Ref. [12], this would be 
k B T < 0.2Ej ~ 2K. A temperature of 0.2-0.4 K should 
be sufficient, and is readily attainable with current cryo- 
genics. 

To summarize, we have demonstrated that the three- 
junction persistent-current qubit can be placed in a 
regime such that the states are determined by the in- 
terference of tunneling paths. For certain values of the 
offset charges, this interference is perfectly destructive, 
leading to a vanishing of the tunnel splitting between 
the two states of the qubit for appropriate values of the 
gate charges and the applied magnetic field. This effect 
should be observable experimentally, as long as the sum 
of the offset charges can be controlled experimentally, as 
we have briefly discussed above. It would certainly be of 
interest to observe the cancellation suggested here. 

This work was supported by NSF through grant 
DMR04- 13395. 
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FIG. 1: Schematic of the circuit for the three-j unction qubit, 
after Ref. [12]. There are two superconducting islands, de- 
noted 1 and 2, whose voltages are Vi and Vi. The three 
junctions in the circuit are indicated by crosses; the i th junc- 
tion has capacitance d and Josephson coupling energy Eji. 
An external flux $ = /$o passes through the circuit. The 
superconducting islands 1 and 2 are also connected to ap- 
plied voltages Va and Vb through capacitors C 9 a and C 3 b\ 
the voltage across these capacitors is V 3 a = Va — Vi and 

VgB=V B -V2. 
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FIG. 2: (Color online) Contour plot of the potential U(fa, fa) 
for the special case a = 1.3 (see text). The horizontal and ver- 
tical axes represent fa and fa. Darker shading means larger 
value of U. For this choice of a=1.3, the state denoted by a 
box with vertical lines in the lower right of the white square 
can tunnel to another state only along the heavy line; for 
other directions, the tunneling barrier is much higher. The 
white square represents the unit cell used in the text. 
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FIG. 3: Contour plot of the splitting 2y/F 2 + \t\ 2 . The hori- 
zontal and vertical axes represent (Qa + Qb)/^ and /. Darker 
shading means larger values of the splitting. Except for the 
value of a, we have used the same parameters as in Ref. [l2l ] 
(see text). The splitting vanishes when {Qa + Qb)/c is an 
odd integer and / = 1/2. 



